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I report on a numerical program, which can be used to calculate any infrared safe three-jet ob- 
servable in electron-positron annihilation to next-to-next-to-leading order in the strong coupling 
constant Os. The results are compared to a recent calculation by another group. Numerical differ- 
ences in three colour factors are discussed and explained. 
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INTRODUCTION 

Jet observables and event shapes in electron-positron 
annihilation can be used to extract the value of the strong 
coupling constant as [ll, l3i ISf • This applies in particu- 
lar to three-jet observables, where the leading-order par- 
ton process is proportional to as- In order to extract 
the numerical value from the LEP data, precise theo- 
retical calculations are necessary, calling for a next-to- 
next-to-leading order (NNLO) calculation. Due to the 
large variety of interesting jet observables it is desirable 
not to perform this calculation for a specific observable, 
but to set up a computer program, which yields predic- 
tions for any infra-red safe observable relevant to the pro- 
cess e^e^ — > 3 jets. Such a task requires the calculation 
of the relevant amplitudes up to two loops, a method 
for the cancellation of infrared divergences and stable 
and efficient Monte Carlo techniques. For the process 
e+e" -^ 2 jets this was done in [J, [s], la, [fl] . In this letter I 
report on a NNLO calculation for three-jet observables in 
electron-positron annihilation. Recently another group 
published results for the NNLO corrections for three-jet 
observables [SJ, |9|, [lOl, [ill. In the calculation presented 
here the methods used are in many parts similar to the 
ones used in [8|, Isl, llO|, [ll | , although I will show that in 
certain points there are important differences. The au- 
thors of [8|, y, |lO|, [ll| made major contributions to the 
development of these methods [5|, [ij, |l3|, [ij, [l5| . 



The numerical results of the two calculations are com- 
pared. The comparison is facilitated by splitting the 
NNLO correction term into individually gauge-invariant 
contributions, such that each contribution is proportional 
to a specific colour factor. For the NNLO corrections 
to e^e~ — > 3 jets there are six different colour factors. 
In three colour factors the two calculations agree {N~'^, 
Nf/Nc, N'i). They disagree in the remaining three colour 
factors (A^c J -^c j ^f^c )■ The numerical differences in 
these colour factors can be traced back to an incomplete 
cancellation of soft-gluon singularities in the calculation 
of refs. [8|, |9|, [lOl, lll|. These singularities require addi- 
tional subtraction terms, which are subtracted from the 
five-parton configuration and added to the four-parton 
configuration. These subtraction terms have a structure 



not present in [lOf and are related to soft gluons. These 
terms occur generically in any NNLO calculation with 
three or more hard coloured partons. 



GENERAL SET-UP 



The perturbative expansion of any infrared-safe ob- 
servable for the process e+e" -^ 3 jets can be written up 
to NNLO as 
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Aa gives the LO result, Bq the NLO correction and Co 
the NNLO correction. The coefficient Co can be decom- 
posed into colour pieces 
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where N^ denotes the number of colours and Nf the num- 
ber of light quark fiavours. In addition, there are singlet 
contributions, which arise from interference terms of am- 
plitudes, where the electro- weak boson couples to two 
different fermion lines. These singlet contributions are 
expected to be numerically small [16|, [iTI, ll8| and ne- 
glected in the present calculation. 

The computation of the NNLO coefficient Co requires 
the knowledge of the amplitudes for the three-parton fi- 
nal state e^e~ -^ qqg up to two-loops 18|, [l9[, the am- 



plitudes of the four-parton final states e~^e —>■ qqgg and 
e'^e'' -^ qqqq up to one-loop [20|, |2l|, [22, 1231 and the five- 



parton final states e''"e~ — > qqggg and e~^e~ —^ qqqqg at 
tree level j24, 25] . Taken separately, the three-, four- and 
five-parton contributions are all individually infrared di- 
vergent. Only the sum of them is finite. However, the 
individual contributions live on different phase spaces, 
which prevents a naive Monte Carlo approach. To ren- 
der the individual contributions finite, several options for 
the cancellation of infrared divergences have been dis- 
cussed, like phase space slicing [20| , sector decomposition 
271 . l28|, a method based on the optical theorem [29[ or 



the subtraction meth od [all 3 . Il3l Il4 lisl Isd . Isil |32| . [SSJ, 
34 ISSl . l36l . ISTI . ISSl . Isg . kd . l4l| . In the present calculation 
I use the subtraction method with antenna subtraction 
terms flSl. 



CANCELLATION OF DIVERGENCES 

To render the individual three-, four- and five-parton 
contributions finite, one adds and subtracts suitable cho- 
sen terms. Schematically, we have 



5 partons : 
4 partons : 
3 partons : 
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Here, dal^ , da^ and dcg are the contributions from 
the original amplitudes with five, four or three final state 
partons. da^^'-' is the NLO subtraction term for four- 
jet observables, containing only three parton tree-level 
antenna functions. At NNLO there are several new sub- 
traction terms required, each of them with a specific 
structure. The term da'^^ contains the four-parton 
tree-level antenna functions. The term da''°°^ contains 
three-parton one-loop antenna functions together with 
tree-level matrix elements and three-parton tree-level an- 
tenna functions together with one-loop matrix elements. 
The remaining terms da'*"''''^'"' , da"''"""*, daP'-°''«=* and 
da'^"^* all contain a product of two three-parton tree- 



level antenna functions. In da 



iterated 



and da 



almost 



one 

antenna function has five-parton kinematics, while the 
other antenna has four-parton kinematics. The former 
subtraction term is an approximation to da^^^, while 
the latter approximates da^ in almost colour-correlated 
double unresolved configurations. In dc/V^oduct \^q\i^ g^j^_ 
tennas have four-parton kinematics. The term da^°f^ will 
be discussed below and is relevant only for the colour fac- 
tors iV2, N° and NfN^. 

The subtraction terms without da*"-''* correspond to 
the subtraction scheme of ref. [10|. For any subtraction 
scheme it is required, that in the three-parton channel the 
explicit divergences cancel, that the four-parton channel 
is integrable over a single unresolved phase space and in 
addition that the explicit divergences cancel and finally 
that in the five-parton channel the integrand is integrable 
over single and double unresolved phase space regions. 
It is easily checked that with the subtraction terms of 
ref. [10| the explicit divergences in the three-parton can- 
cel and I will focus in the following on the four- and 
five-parton channels. 

In the four-parton channel the combination da\ + 
^^NLO jg fj.gg q£ explicit poles. It has been noted 
in ref. [l3| that the combination da'°°P + rfa**'=™*'='' + 



^^product _ ^^aimost jnyolves in the colour factors N^ 
and N^ poles of the form 
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where pii and p2' are the momenta obtained from pi. 
Pi and p2 through a 3 ^- 2 phase space map. ^3 is 
the three-parton tree-level amplitude and X^{l,ij2) a 
three-parton tree-level antenna function. In ref. 10] it 
was claimed that these poles vanish after the azimuthal 
integration over the unresolved phase space. This claim 
is wrong. In the centre-of-mass frame of pi' + p2' with 
Pi' and pi along the positive z-axis, the relevant integral 
is 
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where for x = 2,j we set Cx = cos^j;, s^ = sin^a; and 
6*2 and 9j are the polar angles of partons 2 and j in the 
chosen frame. The integral equals 



In 
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The integral is zero for 9j < 82 but non-zero for 9j > 62- 
In ref. [10] it was claimed that the integral vanishes in 
both cases. As a consequence of the non-zero value for 
9j > 02 the explicit poles do not cancel in the combina- 
tion da'°°P -I- da***^™*'^'' -I- daP™''"'^* - da°'"°"*. The same 
situation occurs also in the colour factor NfNc- 

These poles have a counter-part in the five-parton 
channel. Setting da'^"^* to zero and using a slicing ap- 
proach one observes in the colour factors N^, N^ and 
NfNc a logarithmic dependence on the slicing parameter 
2/rmn = Smin/Q'^- This is showu in fig. [H 

These singularities require an additional subtraction 
term and that is where the present calculation differs 
from the one of ref. [SJ, |9|, [lOl, ll]|. da^°f* is a subtrac- 
tion term related to soft gluons which ensures that the 
poles in the four-parton configuration vanish after inte- 
gration over the azimuthal angle and which renders the 
five-parton configuration independent of ymin ■ The term 
da""-'^* for the four-parton configuration can be taken of 
the form 
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where ^3 is the integrated soft antenna function and pg 
is given by 
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Vcut 


Aa-jet 


B3-jet 


Ca-jet 


0.3 


0.02 


0.13 


-6 ±3 


0.1 


2.12 


34.3 


(2.0 ±0.2) -10^ 


0.03 


7.63 


113.8 


(6.7±0.6)-10'' 


0.01 


15.7 


152.6 


(-1.2 ±0.2) ■ 10^ 


0.003 


27.9 


-6.5 


(-8.1 ±0.5) -10^ 


0.001 


42.4 


-562 


(-21 ±1) ■ 10=* 


0.0003 


61.8 


-1.97 


10^ (-25 ±3) ■ 10=* 


0.0001 


82.9 


-4.36 


10=* (7 ±5) -10=^ 



TABLE I: The LO coefficient A3- jet, the NLO coefficient 
Bz-jet and the NNLO coefficient Cs-jet for the three jet cross 
section with the Durham jet algorithm and various values of 

Vcut- 



I also used the short-hand notation pij2 = Pi + Pi + P2- 
The 6'-function enforces dj > 62 in the specific frame 
introduced above. da^°^^ for the five-parton configura- 
tion is obtained by hfting eq. [5] to the five parton phase 
space. Fig. [T] shows that the sum of all contributions in 
the five-parton channel is now independent of ymin- I 
have checked that in the four-parton channel the explicit 
poles cancel after integration over the unresolved phase 
space. 



NUMERICAL RESULTS 

The numerical program is_build on an existing NLO 
program for e+e^ — > 4 jets [42j- I consider the three-jet 
cross section, where the jets are defined by the Durham 
jet algorithm [43|. The recombination prescription is 
given by the E-scheme. I take the centre of mass en- 
ergy to be vQ^ = rnz- The three-jet cross section is 
expanded as 



Co 
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where ctq is the LO cross section for e^ 



3- jet 



hadrons. 



The coefficients A 



3~jet, ^3— jet 



and C^^jet are given for 



the renormalisation scale /i^ — Q^ and various values of 
the jet defining parameter ycut in table [H The errors 
of Cs-jet are from the Monte Carlo integration. For se- 
lected values of ycut the contribution from the individual 



colour factors to the NNLO coefficient C- 



3~jet 



is shown 



in table lU Finally, fig. [2] shows the scale variation of the 
jet rate defined by 



C^S-jet 
CTtot 

where 



as T 

ZTT 



jet 



(£)^i.3-.,.(g)'c3. 



-jet , 



A3- jet — ^3-jet, 



B3-jet 



B3-jet 



A3-jetAtot, 



C: 



3~jet 



C3- 



3-jet 



B3-jetAtot 



A3-jet [Bt 



■ri-f 



ycut 
0.1 
0.01 
0.001 



ycut 



0.1 

0.01 

0.001 



ycut 
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(1.06 ± 0.02) ■ 10^ (-9.80 ± 0.06) • 10^ 
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-35 ±1 
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(7.09 ± 0.08) • 10^ 
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-0.49 ±0.03 
0.25 ±0.15 
(3.38 ±0.01) ■ 10^ 



Nl 



(1.336 ±0.003) • 10^ 
(1.646 ± 0.002) • 10^ 
(7.41 ±0.01) -10^ 



TABLE IL The contributions from the individual colour fac- 
tors to the NNLO coefficient Cz-jet- 
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44C3MVc + -^±(8C3-ll)iVy 



The renormalisation scale is varied from /z = mz/l to 
/i = 2m z- In this plot the experimental measured values 
are also shown [4^. For values below ycut = 0.001 the re- 
sults of ref. [8|] differ significantly from the ones presented 
here. 



CONCLUSIONS 

In this letter I reported on the NNLO calculation 
for three-jet observable in electron-positron annihilation. 
Particular attention was paid to the cancellation of in- 
frared singularities. I presented numerical results for the 
Durham three-jet cross section. 
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FIG. 1: Dependence of the five-parton contribution on the 
slicing parameter j/mm for the Durham jet cross section with 
= 0.01 in the colour factors iV^^, A''^ and NfNc- "Stan- 



Vcut 

dard" denotes the combination da. 



(0) 



da' 



da' 



^^iterated _ ^^aimost ^ Ug^f^,, ^^^ contribution from da''°fK In 
addition the sum of the two terms is shown. For small values 
of ymin the sum is independent of ymin- 



Durham three-jet rate 




FIG. 2: The scale variation of the three jet rate with the 
Durham jet algorithm at y^^ = mz with as{mz) ~ 0.118. 
The bands give the range for the theoretical prediction ob- 
tained from varying the renormalisation scale from /i = mz/2 
to /i = 2m z- 



